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( Abstract: We study evolution of entanglement of two two- level atoms in the presence of dissi- 

pation caused by spontaneous emission. We find explicit formulas for the amount of entanglement 
as a function of time, in the case of destruction of the initial entanglement and possible creation 
of a transient entanglement between atoms. We also discuss how spontaneous emission influences 
nonlocality of states expressed by violation of Bell - CHSH inequality. It is shown that evolving 
system very quickly becomes local, even if entanglement is still present or produced. 
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The process of spontaneous emission by a system of two - level atoms was extensively studied 
by several authors (see e.g. PEJG3E]). In particular, in the case of spontaneous emission by two 
trapped atoms separated by a distance small compared to the radiation wavelength, where is a 
substantial probability that a photon emitted by one atom will be absorbed by the other, there 
are states of the system in which photon exchange can enhance or diminish spontaneous decay 
rates. The states with enhanced decay rate are called superradiant and analogously states with 
diminished decay rate are called subradiant pQ. It was also shown by Dicke, that the system of 
two coupled two-level atoms can be treated as a single four-level system with modified decay rates. 

Another aspects of the model of the spontaneous emission are studied in the present paper. 
When the compound system of two atoms is in an entangled state, the irreversible process of 
radiative decay usually destroys correlations and the state becomes unentangled. In the model 
studied here, the photon exchange produces correlations between atoms which can partially over- 
come decoherence caused by spontaneous radiation. As a result, some amount of entanglement 
can survive, and moreover there is a possibility that this process can entangle separable states of 
two atoms. The idea that dissipation can create entanglement in physical systems, was recently 
developed in several papers |SJ |S] ■ In particular, the effect of spontaneous emission on de- 
struction and production of entanglement in two - atom system was discussed E3 ID] ■ Possible 
production of robust entanglement for closely separated atoms was shown in Ref. 9 , and the 
existence of transient entanglement induced by this process in a system of two atoms separated by 
an arbitrary distance was also studied [THIHI) . 

In this paper we also concentrate on arbitrarily separated atoms and consider the dynamics of 
entanglement. Similarly as in we take some class of initial states, including interesting pure and 
mixed states, and discuss in details its time evolution as well as the evolution of its entanglement. 
Note that our initial states are different from that considered in Ref. We study also the 

interesting problem how dissipative process of spontaneous emission influences nonlocal properties 
of initial states. Nonlocality of quantum theory manifets by violation of Bell inequalities, and in 
the case of two two-level systems it can be quantified by some numerical parameter ranging from 
for local states to 1 for states maximally violating some Bell inequality. Atomic dynamics studied 
in the paper enables also to consider time evolution of this parameter. 
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The model considered in the present paper consist of two two-level atoms coupled to a common 
thermostat at zero temperature and the reduced dynamics (in the Markovian approximation) is 
given by the semi-group {T t } t >o of completely positive linear mappings acting on density matrices 
(see e.g. |12p. The dynamics takes into account only spontaneous emission and possible photon 
exchange between atoms |13l I14j . and the generator Ljj of {T t } t >o is parametrized in terms of 
the spontaneous emission rate of the single atom 70 and the photon exchange rate 7. In the case 
of atoms separated by an arbitrary distance R, 7 is strictly smaller then 70 and one can check 
that the relaxation process brings all initial states into the unique asymptotic state when both 
the atoms are in their ground states. In contrast to the small separation regime (7 = 70) studied 
in Ref. [5], where the robust entanglement of non-trivial asymptotic states can be analysed, in 
the present case only the transient entanglement of some states can exist. To consider transient 
entanglement we need to know in details time evolution of initial states, not only its asymptotic 
behaviour, so the analysis of possible generation of entanglement is much more involved. 

In this paper we calculate time evolution of an arbitrary initial density matrix. To obtain 
an analytic expression for entanglement as a function of time, we concentrate on the class of 
states which is left invariant during the evolution, and admits explicite formula for the measure 
of entanglement. Next we discuss in details how evolve pure initial states, both unentangled and 
entangled. We show that entanglement as a function of time can behave very differently depending 
on initial conditions: it may monotonically decrease to zero, increase to maximal value and then 
decrease to zero or even it can have local minimum and maximum. In particular, there are states 
for which induced transient entanglement is larger then initial entanglement. Our solution enables 
also to study nonlocality of the evolving system of two atoms. It turns out that the natural measure 
of nonlocality very quickly becomes equal to zero, even if entanglement is increasing in some time 
interval. 
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2. Entanglement and Nonlocality for a Pair of Two-Level Atoms 

2.1. Measure of entanglement. Consider two-level atom A with ground state |0) and excited 
state |1). This quantum system can be described in terms of the Hilbert space Ha = C 2 and 
the algebra 21 a of 2 x 2 complex matrices. If we identify |1) and |0) with vectors (J) and (°) 
respectively, then the raising and lowering operators <r+, er_ defined by 

(1) o-+ = |1><0|, *- = \0){l\ 
can be expressed in terms of Pauli matrices a\ , oi 

(2) (T+ = i (<7i +ia 2 ), a- = - (<Ti - i CT 2 ) 

For a joint system AB of two two-level atoms A and B, the algebra %ab is equal to 4 x 4 complex 
matrices and the Hilbert space Hab — Ha <8> W_b = C 4 . Let £ab be the set of all states of the 
compound system i.e. 

(3) S ab = {p e %ab ■ P > and trp = l} 
The state p G Sab is separable |17j , if it has the form 

(4) p = J2*kp£®Pk, Pi e^A, Pk &£b, A fe >0 and ^ \ k = 1 

fc k 

The set fj^g of all separable states forms a convex subset of Sab- When p is not separable, it is 
called inseparable or entangled. Thus 

( 5 ) £JLb = £ab \ £ab 

As a measure of the amount of entanglement a given state contains we take the entanglement of 
formation ^H] 

(6) E{p) = min ^kE(P k ) 

k 

where the minimum is taken over all possible decompositions 

(7) p = ^Pk 

k 

and 

(8) E{P) = -tr [(tr^P) log 2 (tr A P)] 

In the case of two two-level atoms, E{p) is the function of another useful quantity C(p) called 
concurrence, which also can be taken as a measure of entanglement [151 116j . Now we pass to the 
definition of C(p). Let 

(9) (0* = ((7 2 <8 <7 2 ) P(<72 ® 02 ) 

where p is the complex conjugation of the matrix p. Define also 

(10) p= (^2^1/2)1/2 

Then the concurrence C (p) is given by ^3 E] 

(11) C{p) = max (0,2p max (p) - tip) 

where p m ax(/5) denotes the maximal eigenvalue of p. The value of the number C(p) varies from 
for separable states, to 1 for maximally entangled pure states. 
Consider now the class of density matrices p 



(12) p = 



(° 








\ 





P22 


P23 


P2i 





P32 


P33 


P3i 




Pi2 


PA3 


PAi) 
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where the matrix elements are taken with respect to the basis |1) ® |1), |1) ® |0), |0) (8) |1) and 
|0) ® |0). One can check that for density matrices of the form (12) 

(13) C(p) = |/323 j - \/P22P33 ~ I 1/323 I ~ \J P22pZZ \ 

By positive-definiteness of p, |p23 1 < ^/P22P33, and we have the result: 
Concurrence of a density matrix (12) is given by 

(14) C(p) = 2| P23 | 

As we will show, the class of density matrices given by (12) is invariant with respect to the 
time evolution considered in the paper, and formula (14) can be used to analyse the evolution of 
entanglement of states which initially have the form (12). 

2.2. Violation of Bell inequalities. The contradiction between quantum theory and local real- 
ism expressed by the violation of Bell - CHSH inequality ^5] , can be studied in the case of two - 
qubit system using simple necessary and sufficient condition |2U1 121| . Any state p £ Sab can be 
written as 

(15) / 9=- I 2 «>l2+'r-er<g.II 2 +I 2 <g.s-<x+ 

trim 0~n C8) 0~m I 

\ n,m— 1 / 

where I 2 is the identity matrix in two dimensions, <7i, er 2 , 03 are Pauli matrices, r, s are vectors in 
3 

M 3 and r ■ cr = r j°~j- The coefficients 
j'=i 

(16) t nm = tr (p <j„ ® a m ) 
form a real matrix T p . Define also real symmetric matrix 

(17) U p = TjT p 

where Tj is the transposition of T p . Consider now the family of Bell operators 

(18) B CHSH = a ■ a ® (b + b') ■ a + a' ■ a ® (b - b') ■ a 
where a, a', b, b' are unit vectors in K 3 . Then CHSH inequality reads 

(19) \ty (pBchsh)\ <2 

Violation of inequality (19) by the density matrix (15) and some Bell operator (18) can be checked 
by the following criterion: Let 

(20) m (p) = max (uj + Uk) 

j<k 

and Uj, j = 1, 2, 3 are the eigenvalues of U p . As was shown in |2()l 12 lj 



(21) max tr (p B C hsh) = 2^m(p) 

BcHSH 

Thus (19) is violated by some choice of a, a', b, b' iff m(p) > 1. 

If we consider subclass of the class of states (12) consisting of density matrices of the form 

/0 \ 

(22) p= I P22 P23 ° 

p 32 p 3 3 

\o p 44 / 
then we obtain the following expression for m(p) 

(23) m{p) = max (2 C 2 (p), (1-2 p M ) 2 + C 2 {p) ) 

where C(p) — 2 1 /O23 1 is the concurrence of the state p. Notice that the inequality 

(1-2P44) 2 +C 2 (p) > 1 



is equivalent to 

1/023 1 2 > P44 (1 - P44) 

Let us introduce linear entropy of the state p 

S L (p) = l-tr p 2 

For states (22) 

tr p 2 = p 2 22 + P33 + P44 + 2 \P23 I 2 
so using (p 2 2 + P33 + P44) 2 = 1 we obtain 

Sl(p) = 2 (/O22P33 + P22/O44 + P33/O44 - IP23I 2 ) 

On the other hand 

1/023 | 2 - P44 (/?22 + P33 ) = |y023 1 2 - P44 (1 - P44 ) > 

SO 

P22P33 - 2 ^(p) = IP23| 2 - P44 (P22 + P33 ) > 

and we obtain the following result: 

The states (22) violate some Bell - CHSH inequality if and only if\p23\ > ^75 or P22P33 > 5 Sl(p)- 
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3. Spontaneous Emission and Evolution of Entanglement 

We study time evolution of the system of two two-level atoms separated by a distance R when 
we take into account only the dissipative process of spontaneous emission. The dynamics of such 
system is given by the master equation |131 IT^| 



(24) 



dp 
~dt 



= L D p, p € £al 



with the following generator Lb 



(25) 



L D p= - 



- J2 lki{^-P(T l + -a\a l _p- 

k,l=A,B 



p*W- ) 



where 



(26) 



1 cr±, <r± = ± icr 2 ) 



and 7,4,4 = jbb — Jo, Iab = 1b a = 7 = S7o- Here 70 is the single atom spontaneous emission 
rate, and 7 = (770 is a relaxation constant of photon exchange. In the model, g is the function of 
the distance R between atoms and g — > 1 when R — > 0. In this section we investigate the time 
evolution of the initial density matrix p of the compound system, governed by the semi - group 
{Tt}t>o generated by Lp. In particular, we will study the time development of entanglement of 
p, measured by concurrence. When 7 < 70, the semi-group {T t } t >o is uniquely relaxing, with the 
asymptotic state |0) ® |0). Thus, for any initial state p, the concurrence C{p t ) approaches when 
t — > 00. But still there can be some transient entanglement between atoms [H1E31E]- In this 
section we study in details time evolution of a given initial state p = (pjk)- Direct calculations 
show that the state p(t) at time t has the following matrix elements with respect to the basis 
|1)®|1>, |1>® |0), |0)®|1), |0>® |0) 



Pll 



Pi2(t) =e" 70 * ( P i2 cosh — - pi 3 sinh — ' 

3 + 'yt 'yt 

Pi3(t) =e" 570 * (pi 3 cosh — - pi2 sinh — ' 

p u {t) =e-^ ot p u 



7o " 1 



2^22 - P33) + 



n 2 + il , 1 



2 2 i^ 1L ' 1 



Pu + 7; (P22 + P33)) cosh7* 



9 77 ° 9 P11 - Rep 2 3 sinh 7*) 
7 " T 



7 - T 



lot 



:(P33 - P22) 



7 2 +7o 2 , 1 



P11 + ~x (P22 + P33)) cosh7t 



2 77 ° 9 P11 - Rep23 sinh7i) 

7 - 7 2 



_ e -27o*Jryo^ _ 70t 

7o ~7 2 



- 2(^22 +P33)sinh7t 



7 2 7 2 

(-~ o sinh7t + cosh7t) pn + ilmp 23 + Rep 23 cosh7i 

7o~7 2 



p 24 {t) = [e"5T°* (sinh - 7 cosh |i) + e - ^* (7 cosh |i - sinh p i2 + 
[ e -iTo* (7 sinh |i - cosh |i) + e - ^* (cosh |i - 7 sinh p 13 + 
e~~* (P24 cosh -t - p 34 sinh -i) 



p 34 (t) = [e"5T°* (7 sinh - cosh ^i) + e 2 * (cosh - 7 sinh ^i)] p i2 + 

[ e -f 7o* ( sinh l t _ 7 cosh 2 t ) + e -^t ( 7 cosh 2 t _ sintl pi3+ 



e 2 * (p34 cosh -t - p 24 sinh -t) 



7o ~7 2 



7o -7 

(P22 + P33) cosh7i + 2Rcp 23 sinh7i 



' 4 77o . , 2(7 2 +7^) N 
( 2 _ " 2 sinh 7 i - w 2 _ 2 ° ; cosh 7t) p n 



7o -7" 



The remaining matrix elements can be obtained by the Hcrmiticity condition. One can simply check 
that the classes of states (12) and (22) are invariant with respect to the above time evolution. In 
particular, in both cases 



P23(t) = e 7ot [Rep 23 cosh7i + zlmp 23 - - (p 22 + p 33 ) sinl^t] 



C(p(t)) = 2e 7o * I (Rcp 23 cosh7i + iIm / 9 23 - ^ (p 22 + p 33 ) sinh 7* ) | 



(27) 

So we obtain the result: 

For initial states (12) or (22) the concurrence at time t is given by 
(28) 

where pjk are matrix elements of the initial state. 

Let us consider pure initial states ^ € C 4 belonging to the class (12). The most general pure state 
of this type can be written as 

(29) * = cos^cos^ll) ® |0) + sin cost/; e ie |0) <g> |1) +sin^e lH |0) <g> |0) 
with 0, ip e [0, f ], 6, He [0, 2n]. Using (14) we see that concurrence of (29) is given by 

(30) C(^) = cos 2 ip sin 20 

By (28), the time evolution of this initial concurrence is described by the following function 

(31) C(p{t)) = e~ 7ot cos 2 ip I sin 20 cos0 cosh7i — sinh7i — i sin 20 sin 6 | 
The function (31) is simple to analyse when = or ^. In this case 

(32) * = e 40 cosV>|0)<g)|l) + e e sin V> |0) <g> |0) 
and 

(33) C{p{t)) = e~~< at cos 2 V sinh 7 t 
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From the formula (33) we see that initial concurrence equal to zero, increases in the time interval 
[0,imax], where 

(oa\ + 1 7o+7 
(34) t max = — In 

2 7 70-7 

to the maximal value 

2 ., 7 ( 7o + 7 



(35) C max = cos tp 

7o-7 V7o - 7 

and then asymptotically goes to zero (see Fig. 1) 




Fig. 1. C(p(t)) for initial states (32) with tp = and different values of 7/70 

If is arbitrary, we can put for simplicity -0 = 0. Then 

(36) * = cos0|l) <g> |0) + sin0e ie |0) <g> |1) 
and C(\I/) = sin 20. The evolution of this initial concurrence is given by 

(37) C{p{t) ) = e I sin 20 cos 6 cosh jt - sinh jt - i sin 20 sin 6 | 

Depending on values of and O, the function (37) can be strictly decreasing to zero, or can have 
one maximal value for some t, or even can have one minimal and one maximal value. Let us discuss 
all these possibilities by choosing some special initial states from the class (36). 

a. Let 9 = 0. Then 

(38) * = cos0|l) ® |0) +sin0|O) ® |1) 
and 

(39) C{p{t)) = e~ 70t I sin 20 cosh jt - sinh jt I 
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The function (39) is decreasing to zero in the interval [0, t m i n ] where 

1 , 1 + sin 20 

40 tmin = ^ In- 

2 7 1 — sin 20 

Then in the interval [t m in, tmax], with 

(41) t 1 ln (j±£^fe±l) 

( j W -2 7 in (l-sin20)( 7o - 7 ) 

(39) increases to the maximal value 

, . _ 7 |cos20| f (l + sin20)( 7o + 7 )V^ 



"0 



^2 V(l-sin20)( 7o -7) 

For i > i max , (39) goes asymptotically to 0. For sufficiently small initial concurrence C max > C(*f>) 
but for larger entanglement of the initial state, the maximal entanglement produced during the 
evolution is smaller then C(\I/) (see Fig. 2. below). 



C( P {t)) 
0.5 



. 4 



0.3 



0.2 



0.1', 




7 8 7o* 



Fig. 2. C(p(t)) for initial states (38), with = tt/40, tt/20 and 7/70 = 0.75 



b. Let 6 = 7T. Then 

(43) 

and 

(44) 



* = cos0|l) <g> |0) - sin0|O) <g> |1) 



C(p(t)) = e" 70 * sin 20 cosh 7 i + sinh 7 i 



If sin 20 > ^ then function (44) is monotonically decreasing to 0. On the other hand, if sin 20 < ^ 
then at time 

(45) t ~ 1 ln ( 1 - sin2 ^)(7Q+7) 

1 j maX_ 2 7 ln (l+sin20)( 7o - 7 ) 
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(44) attains local maximum 

_ 20. 

7 I cos 20| /(l^sin20)(7o+7)\ ^ 
( j v^o^ l(l + sin20)( 7o - 7 )j 

Cmax is always greater then initial concurrence C(^f) and 

(Af\ n 1 + sin 20 

(47) C max -> when 7 -> 70 

and 

7 7 

(48) C max — ► — 10/ien sin 20 — > — 

7o 7o 
Thus, for cntanged pure initial states (43) the dissipative process of spontaneous emission increases 
entanglement, provided that the initial entanglement was smaller then ^ (see Fig. 3). 




. 4 



7 8 7o* 



Fig. 3. C(p(t)) for initial states (43) with C(*) = 0.1, 0.4, 0.7, 1 and 7/70 = 0.75 



* = cos0|l) <g> |0) +i sin0|O) <g> |1) 



c. Let 9 = tt/2. Then 
(49) 
and 

(50) C{p{t)) = e-^ at j sinh jt + i sin 20 | 

One can show that if | sin40| < 7/70 then (50) achieves local minimum at 

Or ~2 ~2 



(51) imin 

and local maximum at 



I 70 cos 40 - y 7 2 - 7^ sin 40 



2 7 



In 



7o-7 



(52) 



I 70 cos 40 + y 7 2 — 7q sin 40 

^max = 7: In 

27 70-7 



11 



with the corresponding values of concurrence 

70 cos 40 — J 7 2 — 7g sin 2 40 



(53) C„ 
and 

(54) C a 



20 

2". 



70-7 



1 



7 2 cos 40 — 7 ^/7 2 — 7q sin 2 - 
2(7o 2 -7 2 ) 



7o cos 



40+ ^7 2 - 7o sin2 ' 
7o - 7 



7 2 cos 



7 2 - 7 2 sin 2 40 



2(7 2 -7 2 ) 



For other cases, the function (50) is monotonically decreasing to (see Fig. 4. below). 



C(p(t)) 




1 2 " 3 4 " " 5 6 7 8 lot 

Fig. 4. C(p(t)) for initial states (49) with = tt/20, tt/40 and 7/70 = 0.75 
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4. Evolution of Nonlocality 

Nonlocality of quantum theory manifesting by violation of Bell inequalities is strictly connected 
with the existence of entangled states. It is known that every pure entangled state violates some 
Bell inequality (see e.g. [221) - But for mixed entangled states it is no longer true ^7j. So it is 
interesting to discuss how dissipative process of spontaneous emission influences nonlocal properties 
of initial states. For simplicity we restrict the class of states considered below to density matrices 
of the form (22). For that class we can apply the results of Sect. 2.2. 

As the initial states we take the states (43). Note that at time t, the state pit) will have the form 
(22). The initial entanglement is non-zero, so violation of some Bell - CHSH inequality occurs at 
time t = 0. What happens during the evolution? Consider the inequality 

(55) P2 2 (t)p 33 (t) > \s L (p{t)) 

which is sufficient to nonlocality of the state pit). Observe that Sl(/9(0)) = and S^pit)) is 
increasing in some time - interval. On the other hand, 

P22(0)p 33 (0) > |p 23 (0)| 2 > 

and P22{t)p3z{t) asymptotically goes to 0, so there is some non - empty interval < t < t\ for 
which the inequality (55) is satisfied. Thus for < t < ti, all states pit) will still have nonlocal 
properties. We may also introduce the time t n after which nonlocality is lost. To this end, besides 
ti consider t2 such that 

lP23(t2)\ = ^ 

Then 

(56) t„ = max (ti.t 2 ) 

By the results of Sect. 2.2, all states pit) for t > t n will admit local hidden variable model. To 
illustrate this concept, consider as initial states 4 ,+ and (symmetric and antsymmetric states) 

* ± = 7f (|1>®|0>±|0>®|1» 
Let p^it) denote corresponding density matrices at time t. Then 



\P2s(t)\ = Je"^)' 



and 

We see that 
Thus 



pUt)pUt) = \ e- 2 ^^\ S L { P ±{t)) = 2 (e-(70±7)* _ e -2(70±7)*) 



15 1 In 2 

In ■ 



7o ± 7 4 70 ± 7 

1 In 2 

tn 



7o ± 7 2 

Note that for antisymmetric state ^~ and 7 close to 70, t n goes to infinity. 

It is also interesting to study in more details time evolution of the measure of nonlocality which 
may be defined as follows 

(57) nip) — max ( 0, m(p) — 1 ) 

As is well known, m(p) < 2, so < nip) < 1 and larger value of nip) gretar then 1 means violation 
of CHSH inequality to a larger extent. Since m(p) = 2 for maximally entangled pure states which 
maximally violate CHSH inequalities, for them n(p) = 1. To obtain analytic expression for n(p(i)) 
we can utilize formula (23). Formula (23) is further simplified if we take such initial states that 

(58) \p2s(t)\ < ^= 
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for all t. This condition can be achieved if 



~< J r 

7o V2 



Then inequality (58) is satisfied and m(p(tj) > 1 if and only if 
(59) (l-2p u (t)) 2 + C 2 (p(t)) > 1 



and to study the time evolution of nonlocality, we only need to know time-dependence of the left 
hand side of (59). If we take initial state with sin 20 < ^ < -^=, then during the time evolution 

C(p(t)) increases to the maximal value C max < -±= in the interval [0,i max ]. At the same time, 



(60) (1 - 2p 44 f = (2e" 7ot (cosh 7 i + sin 20 sinh 7 t) - 1 



decreases so fast that left hand side of (59) is a decreasing function of t (sec Fig. 5, 6 below). 
Thus we obtain the result: 

For initial states (43) and sin 20 < ^ < nonlocality of p(t) given by n(p(t)) decreases during 
the time evolution even if the entanglement increases. 
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n(p(tj) 

1 — - 

0.8 
0.6 
0.4 
0.2 




0.02 0.04 0.06 0.08 0.1 Tot 

Fig. 6. n(p(t)) for initial state (43) with C(*) = 0.3 and 7/70 = 0.7 

Numerical analysis indicates that for other initial states (43) for which m(p(tj) is defined by the 
whole expression (23), the nonlocality n(p(tj) also monotonically goes to irrespective of the 
evolution of entanglement (see Fig. 7 below). 



n(p(t)) 




0.25 0.5 0.75 1 1.25 1.5 1.75 2 lot 
Fig. 7. n(p(t)) for initial states (43) with C(*) = 0.8, 0.9, 1.0 and 7/70 = 0.7 
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